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Supporting human life on a lunar base uill require grounng many different food crcps. This paper 
investigates the gyowth dynamics of four crops (wheat, soybeans, potatoes, and lettuce) for general 
similarities and differences, along uith associated material flows of the gases, liquids, and solids in 
a lunar farm. The human dietary requirements are compared uith the protein, carbohydrate, and Upid 
contents of these hydroponically grown, highproductiiity cnps to derive a lunar farm diet. A simple 
and general analytical model is used to calculate the mass fluxes of CO 2 , H 2 O, HNO:^ and O 2 during 
the life cycle of each of the four crops. The resulting farm cnp areas and corresponding biomass 
production rates are given. One significant conclusion of this study is that there is a “lipid problem’' 
associated uith the incorporation of these four crops into a viable diet. 


INTRODUCTION 

Following the return of our astronauts to the lunar surface 
around the turn of the twenty-first century, an outpost for 
temporary habitation could evolve into a permanently occupied 
base on the Moon {Ride, 1987). The major human life support 
needs will have to be met at increasingly self-sufficient rates 
during this evolution. The pathways leading to a lunar farm are 
yet to be defined in the habitat development scenarios. 

Human diets for a lunar base can be provided with hundreds 
of foods. Here, however, we will focus on four crops studied in 
the NASA Controlled Ecological life Support Systems (CELSS) 
Program: lettuce, potatoes, soybeans, and wheat. Substantial data 
have been generated on the response of these crops to variables 
important in future space agriculture such as near-maximally 
achievable planting density, light intensities and schedules, and 
atmospheric CO 2 levels. Additional experimental data for these 
crops were received in 1987 through personal communication 
with CEUSS researchers B. Bugbee, C. Mitchell, D. Raper, 
R. Wheeler, and S. Schwartzkopf Information received included 
environmental conditions for both the aerial and root plant parts 
in particular high-yield experiments. Figure la shows the com- 
position of the edible portions of lettuce, potatoes, soybeans, and 
wheat in terms of the three major food types, protein, carboh- 
ydrate, and lipid. 

To incorporate these crops into a farm, we consider the dietary 
needs that must be met by the candidate crops. Figure lb shows 
the protein, carbohydrate, and lipid requirements of two standard 
satisfactory diets. More detailed dietary breakdowns, such as 
essential amino acids, fatty acids, and vitamins are beyond the 
scope of this study. Even though each diet provides 2700 kcal per 
day per person, the relative fractions of calories obtained from 
proteins and lipids are different. By comparing the compositions 
of the crops (Fig. la) with those of the diets (Fig. lb), a lipid 
problem becomes evident. 


The lipid problem arises because both standard diets contain 
more lipid than protein. Diets with lower lipid than those used 
here might be desirable {Roberts, 1988). Because none of the 
four crops contains more lipid than protein, any allotment we 
make using these crops to fulfill the total lipid requirements will 
concomitantly have an excess of protein. Waste such as this would 
be detrimental to a space agriculture prescribed by energy and 
mass constraints. 

CROP MODEL DEVELOPMENT 

Simulation models help us conceptualize and design new 
systems by using a mathematical framework to assemble 
components for investigating specific system-level issues. Previous 
work along these lines developed a model (called BLSS) for a 
CELSS that grows wheat as the sole crop {Volk and Rummel, 
1987; Rummel and Volk, 1987). BLSS can be used to track the 
flow of carbon, hydrogen, oxygen, and nitrogen through the 
various processes in a CELSS because it contains the stoichiome- 
tries for various compounds such as plant protein and human 
urine. The model grows wheat in a variety of planting schemes, 
with different numbers and sizes of simultaneous batches. 
Different schemes produce different magnitudes of fluctuations in 
the standing biomass and in the buffer mass reservoirs of C 02 > 
H 2 O, HNO 3 , and O 2 . 

Here we extend this approach to include lettuce, potatoes, and 
soybeans also. Figure 2, along with the model results still to be 
discussed, shows selected and typical data for the growth of the 
edible and inedible parts (to humans) of each crop. A breakdown 
of biomass into edible and inedible parts is fundamental in a 
CEI.SS because of the consequent separation of material flows. 

Many crop growth curves prominently show an S-shaped or 
sigmoidal curve typical of biological systems. The logistic 
differential equation dC/dt = rC(l-C/K) imitates this S-shape of 
exponential growth followed by a leteling off. The term C is 



526 2nd Conference on Lunar Bases and Space Activities 


biomass, t is time, r is growth rate for the purely exponential part 
of the system, and K is a “negative feedback” from the growth 
process itself, an environmentally modifiable but inherent (genet- 
ically based) slowing of the total growth rate (dC/dt) by the 
approach of the crop to its mature size. The logistic equation thus 
contains some biologically meaningful parameters and is chosen 
to represent the growth of the inedible crop parts. 

The equation for the edible crop parts must be somewhat 
differently structured. The edible cells, like the inedible ones, 
reproduce, so the total edible growth is set proportional to the 


edible mass. Furthermore, the nonphotosynthesizing edible parts 
(except for lettuce; see below) grow using products from 
photosynthesis by the inedible parts (the leaf mass); therefore, 
the inedible biomass (Mir^.^) should also appear in the edible 
equation. Also, the edible growth occurs substantially after the 
beginning of the inedible growth (see Fig. 2), so a switch-on time 
(t*) is used in the formulation for edible growth. The edible 
biomass (M^-^j) is assumed to be equal to zero before t* and to 
start its growth at t* with minimum edible mass (E^in). With 
these considerations we write 
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The parameters t and t* are in units of time, while and 
in tirne"^ and the remainder in mass (see Table 1). For wheat, 
soybean, and potato we use equations (la) to (Ic). Because the 
edible and inedible parts develop together, the parameter t* is 
defined differently for lettuce. Mitchell et al. (1986) found that 
the growth rate increases by more than a factor of two at about 
1 1 days; therefore we define r^j 2 and r^d,i for t > t* and t < t*, 
respectively. The equations become for lettuce 
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Fig. 1. (a) Compositions of lettuce leaves, potato tubers, soybeans, and 

wheat berries for typical high-yield hydroponic growth experiments. Data 
provided by CEI.SS researchers C. Mitchell (lettuce), R. Wheeler 
(potatoes), D. Raper (soybeans), and from Bugbee and Salisbury ( 1988, 
wheat). The balancing components of fiber and ash are not shown, 
(b) Compositions of two possible diets. Diet A is from the 1980 
Recommended Dietary Allowances and Estimated Safe and Adequate Daily 
Dietary Intake, using American Heart Association recommendatioas of 35% 
of food kcal from fet {Krause and Mahan, 1980). Diet B uses the NIH 
recommendations (C. Mitchell, personal communication, 1988) of 0.5 g 
protein per day per lb of body mass and using lower value of the 
recommended 30-50% of nonprotein food kcal as lipid to give low er lipid, 
higher protein diet to contrast with diet A. Both diets are approximately 
for a 1 55*lb individual having 2700 kcal per day. 


These models were run in a computer program and the results 
generated were compared to the experimental crop data. 
Adjustments were made to the parameters until the models agreed 
reasonably with the data. The parameters used for each crop are 
listed in Table 1, while the model outputs are shown in Fig. 2. 

The output curves demonstrate that it is relatively easy to 
represent the data with a model whose parameters have some 
fundamental biological meaning. Table 1 lists the actual planting 
mass for the crops, but we need to investigate further the data 
at t = 0 to determine whether they corre^ond to the initiation 
of the crop from seed or tissue or to the transplanting time after 
initial seeding growth. Some further adjustment might be 
necessary to account for the physical meaning of time t = 0. 

Additional refinements are possible. Better fits to the growth 
curves shown for wheat and potatoes in Fig. 2 are obtainable. 
More importantly, the model parameters, such as growth rates 
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Flc. 2. Models of crop growth using parameters from Table 1, compared to crop growth data, (a) Lettuce data are from Mitchell et at. ( 19«>) 
at 1000 ppm CO 2 and 450 jumol/m^-sec of PPF (data were given per plant and adjusted here to yield leaf production of 60g/m d (C. Mitchell, 
personal communic-ation 1987)). (b) Potato data are from Wheeler and Tthhitts (1987) for dry mass produrtion under 24-hour continuous light 
at 300umol/m-=-sec PPF (assume 5 plants per m^). (c) Soybean data are from D. Raper (personal communication. 1987) grown at 700 Mmol/m - 
sec PPF and 400 ppm CO, (data were interpolated by D Raper to be in equal time intervals), (d) Wheat data are from B. Bugbee 
communication, 1987) for plants grown at 1200 Mmol/m^ sec and 1200 ppm CO^ (sec also Bughee and Salisbury, 1988). l>ata represent inthvtdual 
erowth experiments, not neces,sarily the maximum yields ever obtained. Model parameters were not adjusted to achieve exact fits to growth data, 
rather to demonstrate the utility of equatioas (1) and (2) in providing a relatively simple method of generating growth curve-s to determine gas 
and fluid fluxes applicable for including plants in systems m(xlels. 


TABUi 1 . Parameters for crop models 


Parameter 

Wheat 

Soybean 

Potato 

I>etluce 

rinc-dCtlay ') 

0-09 

0.10 

0.06 

same as 

rrd(<lay ') 

0.17 

0.10 

0.30 

r,.d.i =0.2, r^.j,2 = 0.5 


3700.0 

1300.0 

lOOO.O 

1000.0 


2500.0 

1100.0 

4000.0 

5000.0 

Emin 

80.0 

80.0 

400.0 

X 

^ined.o 

150.0 

20.0 

25,0 

X 

Mcu,<i 

0.0 

0.0 

0,0 

2.0 

todays) 

45.0 

45.0 

40.0 

11.0 


Units for K,nc-ti> Emm, Mintxi..., are g dry mass m \ 


(fjS) and ultimate biomass (K^sS), are not constant, but are 
functions of environmentaJ conditions. A reasonable approach 
could be to develop these parameters along the lines of classical 
mathematical treatments of photosynthesis, such as in Gates 
(1980), wherever possible. That way the data would not be used 
for fitting, but rather for model validation. Transpiration sub- 
models and the relationships between atmospheric pC 02 , 
humidity, nutrient uptake, and biomass growth need to be 
developed for investigation of the various design tradeoffs 
between energy, mass, and volume. The models shown here 
would serve as a basis for further developments. 

Volk and Rummell (1987) listed formulas for protein, 
carbohydrate, lipid, fiber, and lignin that can be placed into 


a-3 
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balances equations containing carbon, hydrogen, oxygen, and 
nitrogen. It is therefore possible to calculate the uptake of CO2, 
H2O, and HNO3, and the production of O2 by the crops. These 
compounds vary as a function of fractional distribution of protein, 
carbohydrate, lipid, fiber, and lignin in the biomass. Table 2 shows 
the mass balances for the four crop models. For example, note 
the substantial differences between soybean and wheat in the CO2 
required and the O2 produced per gram of edible biomass 
produced. This difference is due primarily to the difference in 
lipid content. There are corresponding differences in the fluxes 
of these materials between the crops and their environments. 
These fluxes are important in the design of engineered hardware 
for the various crops. 

The balances in Table 2 were used with the crop growth 
models to calculate the fluxes of CO2, H2O, HNO3, and O2 during 
growth; these fluxes are shown in Fig. 3 . Note the different curves 
for the crops. Such curves will be produced during the actual 
operation of a CELSS (e.g., if CO2 will be monitored and 
maintained at desired levels in the crop’s atmosphere, the amount 
of CO2 injected to maintain these levels will be known). Due to 
the characteristic patterns of these fluxes, it is possible to relate 
this information to the monitoring system for the state of the 
whole crop. Note that these curves assume a constant percentage 
of protein, carbohydrate, lipid, fiber, and lignin for the edible and 


TABLE 2. Mass balances for crop models. 


Mass Types 

Wheat 

Soybean 

Potato D'ttuce 

Edible Mass fractions 





Protein 

0.17 

045 

0.13 

0.26 

Digestible Carbohydrate 

0.78 

0.30 

0.84 

0.12 

Upid 

0.02 

0.25 

0.00 

0.06 

Fiber 

0.03 

• 

0.03 

0.56 

Lignin 

0.00 

« 

0.00 

0.00 

Fluxes During Edible Biomass Production 

( g per g dry biomass) 


CO2 (in) 

1.62 

2.10 

1.57 

1.82 

H20(in) 

0.59 

0.66 

0.58 

0.57 

HNO3 (in) 

0.13 

0..54 

0.10 

0.20 

O2 (out) 

134 

2.11 

1.26 

1.59 

Inedible Mass Fractions 





Protein 

0.09 

0.17 

0.19 

0.11^ 

Digestible Carbohydrate 

0.14 

0.80 

0.30 

0.11^ 

Lipid 

0.00 

0.03 

0.00 

0.00^ 

Fiber 

0.72 

* 

0.45 

0.78^ 

Lignin 

0.05 

* 

0.06 

0.00^ 

Fluxes During Inedible Biomass Production (g perg dry biomass ) 


C02(in) 

1.72 

1.63 

1.75 

1.68 

H20(in) 

0.56 

0.59 

0.56 

0.55 

HNO^Cin) 

0.07 

0.13 

0.14 

0.08 

O2 (out) 

L35 

1.36 

1.45 

1.32 


Fiber and lignin were included in the soybean carbohydrate data. 
^ Values assumed by T Volk. 


1 

CO 2 flux 

2 

O2 flux 

3 

Dry blomaai flux 

4 

H2O flux 

5 

HNO3 flux 






Fig. 3. Fluxes of CO2, metaboUc H2O, nutrient HNO3, O2 produced, and total dr>^ weight biomass (edible plus inedible) for the four crop mcxlels: 
(a) lettuce; (b) potato; (c) soybean; and (d) wheal. Note different units for the different crops. Fluxes are from the models of Fig. 2 using the 
stoichiometries of Table 2. 
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inedible during their respective growths. That this is clearly not 
the case is seen in the decrease in leaf N during the seed growth 
in the hydroponic wheat (^Bugbee and Salisbury, 1988). A next 
step here would be to let this N change represent a decrease in 
the protein of the edible parts during the late state of growth 
and to see how much this decrease affects the CO 2 , H 2 O, HNO 3 , 
and O 2 fluxes. 

LUNAR FARM DISCUSSION 

The crops can be incorporated into a collective model for the 
entire farm, assuming the relative areas and volumes for each crop 
are known. We now assemble the four crops into a diet follov^ng 
a particular logic. We first assume that a person could consume 
10 g of dry biomass of lettuce leaf per day. Furthermore, to use 
all four crops and take advant^e of the complete protein created 
by the combination of grains (wheat) and legumes (soybean), we 
assume equal contributions from potatoes, wheat, and soybean to 
meet the daily protein requirements. After satisfying the protein 
requirements, the next critical component is lipid. The only crop 
with substantial lipid is soybean, so additional soybean is added 
to bring the total lipid up to the target values for the two diets. 
All these results are summarized in Table 3 


The protein and lipid requirements are now satisfied, but 
carbohydrate is still short. Potatoes have a significant fraction of 
carbohydrate, with a ratio of carbohydrate to protein approxi- 
mately the value required by the diets. The final step in forming 
the diet, therefore, is to add potatoes until the target value for 
carbohydrate is reached; but this adds stiU more protein. As seen 
in Table 3, the mix of crops to yield 100% of the target values 
for protein and lipid results in an excess of protein, with total 
protein now about 400% and 250% of the respective requirements 
for diets A and B. 

By considering the areas required to grow each crop, the total 
farm area for the life support system can be estimated ( see 
Table 4). The per-area productivity for each crop used in this 
computation was taken from the data used in Fig. 2. Note that 
s(^me of these crops have been grown at higher productivities; 
wheat, for example, has been grown at double the productivity 
shown by increasing the light level {Bugbee and Salisbury, 1988). 
Thus higher light levels might yield still higher productivities, 
light will probably be a useful control parameter for temporarily 
increasing the yields following crop failure or equipment 
downtime when storage reservoirs need increased rates of 
replenishment. Thus the productivities shown in Table 4 were 
deliberately chosen not to be the maxima. For one thing, the 


TABLE 3 . Assembl>^ of a lunar farm diet with four crops. 


Diet A Dirt B 


Crop 

Protein 

Carbo- 

hydrate 

Lipid 

Dry 

Mass 

Protein 

Carbo- 

hydrate 

Lipid 

Dry 

Mass 

Lettuce 

2.1 

1.0 

0.5 

10.0 

2.1 

1.0 

0.5 

10.0 

Potato 

18.7 

118.0 

0.0 

143.8 

25.8 

162.8 

0.0 

198.5 

Soybean 

18.7 

11.0 

10.6 

42.5 

25.8 

15.2 

14.7 

58.6 

Wheat 

18.7 

87.7 

2.2 

114.7 

25.8 

121.0 

31 

158.3 

Soybean 

160.9 

94.9 

91.7 

365.7 

107.7 

63.6 

61.4 

244.8 

Potato 

11.1 

70.4 

0.0 

85.4 

8.6 

54.4 

0.0 

73-8 

Total 
% taiget 

230.2 

4ll 

383.0 

100 

105.0 

100 

762.1 

195.8 

253 

418.0 

100 

79.7 

100 

744.0 


*Tareet values for protein are S6 g day ' for diet A and 77.S g day for diet B i ) 

' Target values for lipid are 1 05 g day ’ for diet A and 79.7 g da>’ for diet B ( see Fig. I ). 

* Target v-aliics for carbohydrate are g day ' for diet A and 418 g day for diet B (see Fig. ). 
All values except percentages are in g person * day'V 


Rationale 


Assume 1 0 g dry mass 
person' ‘ day ’ 

Assume 1/3 taiget 
protein* supplied 
Assume 1/3 taiget protein 
supplied 

Assume 1/3 taiget protein 
supplied 
Add soy 

until lipid target^ 

Add fK>tato until 

caibohydrate target^ 


TABIT 4. Illustrative crop areas for the lunar farm 





Diet A 



Diet B 


Crop 

Prtxluctivit)^ 
of edible 
mass 

g day ' 

Required 
edible 
pnxluction 
gperson ’ day 

(Growing 
area per 
person 

Cirowing 
area for 1 2 
people 
m^ 

Required 
edible 
prixluction 
g person ’day 

Grt)wing 
area per 
person 

Cirowing 
area for 1 2 
people 

m'^ 

U*ttuce 

Potato 

So>bean 

Wheat 

Total 

60 

27 

11 

30 

10 

229.2 

408.2 
1 14 7 
762.1 

0.2 

8.5 

37.1 

3.8 

49.6 

2.4 

102,0 

445.2 
45.6 

595.2 

10 

272.3 

303.4 
158.3 
744.0 

0.2 

10.1 

276 

5.3 

432 

24 

121.2 

331.2 

636 

518.4 


•pr.Hlucltvi.ics arc illusirativc only, no. maximum for each crop WTica., for example, hxs been powm a-s^hij^ as «)g m ■ 
day ' bu. the value of id is used here so higher illumiiution could be ased as a eonlntl lo allow for hii^icr prixJuaion 
unk-r unusual eireumslances I. will be a.s.sumed .ha. dK- .Hher crops are similar in having higher pnxluetm.ies in conditions 
Still to be investigated 

^Noie this amount of sovticans creates a wasteful excess of edible protein ( see Fable 5 ) 
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maxima are not yet known. For another, the production rates 
during normal operations will be less that the maxima to allow 
the system to be controlled when storage reservoirs need to be 
readjusted. The productivities used here are representative of 
hydroponic crop yields that could be accomplished with today’s 
technology. 

As apparent in Table 4, using all the preceding calculations with 
attendant assumptions, most of the area of a lunar farm will be 
dedicated to soybeans (75% for diet A, 64% for diet B). This is 
a direct result of using soybeans to match the lipid requirements. 

CONCLUSIONS 

We have shown that a simple, generic crop model can represent 
the growth of four different candidate crops for Controlled 
Ecological Life Support Systems, providing mass fluxes associated 
with growth for any whole-system CELSS model. An initial 
simplicity is desirable because the model will tend to quickly 
become more complex when it incorporates additional refine- 
ments, particnilarly sensitivities to environmental variables. There 
is every reason to expect that a generic model like the one 
demonstrated here will be useful in constructing a new model 
system for studying the dynamics of a ^ace farm. 

An important problem exists in attempting to combine the four 
crops of lettuce, potatoes, soybeans, and wheat into an adequate 
diet. Besides being bland, there will be a serious overproduction 
of protein. Either diets with much lower lipid content than those 
shown must be designed and approved, or other crops with a 
higher lipid-to-protein ratio should be included. Rapeseed, for 
example, is about 50% lipid and about 20% protein; peanuts can 


be grown with as high as 54% lipid and as low as 21% protein 
(C Mitchell, personal communication, 1988). If these crops were 
used to satisfy the lipid requirements, protein excess could be 
avoided. Unfortunately, little is known about the behavior of these 
crops in high production hydroponics. We recommend systematic 
crop growth experiments aimed at a balanced diet with minimal 
waste. 
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